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Abstract 

We present the hamiltonian study of super Yang-Mills quantum mechanics (SYMQM). The 
recently introduced method based on Fock space representation allows to analyze SYMQM 
numerically. The detailed analysis for SYMQM in two dimensions for SU(3) group is given. 

1 Introduction 

Supersymmetric Yang- Mills quantum mechanics are very interesting models since they emerge 
in different areas of physics. The general, not necessarily gauge, supersymmetric quantum 
mechanics have been studied first as a laboratory of supersymmctry [1] where in particular the 
exact solution for D — 2, SU(2) case was given. By definition SYMQM are AT — 1 super Yang- 
Mills field quantum theories reduced from D = d+ ltoD = + l dimensions. Supersymmetry 
requires the space-time dimension to be D = 2, 4, 6, 10 with M = 2, 4, 8, 16 supercharges in the 
resulting quantum mechanics respectively. The rotational symmetry and gauge invariance of 
the original theory become now the internal Spin(d) and global SU (N) symmetry. The physical 
states become now the SU(N) singlets. We denote the spatial components of gauge field A l a (t) 
by x\ and their conjugate momenta by p l a , \x x a ,p^ ~ ^ lJ ^ab- The hamiltonian is then [1] 

H = \vlV l a + -/{faboXlxif + Hp , (1) 



where H F = -^gfabc^a x b^a)3^c f° r D = 2,10 and ?9 are real spinors obeying } = S al3 S a b, 

a, (3 — l,...,J\f or Hf = igf a bc^2 x l^a^c f° r D = 4,6 and ■& are complex spinors obeying 

— Tho . , 

algebra {T l ,T^} = 28^. 



= S al3 S ab for a, /3 = 1, . . . , =f . The are matrix representation of an SO{d) Clifford 
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The growing interest in these models is due to the BFSS ( Banks, Fischler, Shenker, Susskind 
) conjecture [4] where the N — > oo limit of Eg.(l) is argued to describes M-theory in the infinite 
momentum frame. This stimulated further work on asymptotic form of the ground state of 
D=9+l, SU(2), SYMQM [8] and the analysis of Wittcn index of (1). The index does not vanish 
only in D=10 where it is equal to 1 [9,10,11,12]. Despite the relevance to M-theory SYMQM have 
been studied earlier in different context. The bosonic part of (1) was discovered in pure Yang- 
Mills theory in the zero volume limit [2]. Later on it appeared as a regularization describing 
the quantum supermembrane [3]. The detailed study of the hamiltonian (1) shows that in 
bosonic sector the potential is confining and there is no continuous spectrum [6]. If however 
the supersymmetry is turned on then there are bound states in fermion rich sectors as well as 
scattering ones [7]. 

The only exact solutions of (1) existing in the literature are for D=l+1, SU{2) [1] and its 
generalization for arbitrary SU{N) [5]. Therefore any numerical approach is of interest. 

The plan of this paper is the following. In section 2 we briefly outline the method used to 
study the models just described and quote existing results in D=l+1,3+1,9+1 for SU(2) group. 
In section 3 and 4 we study general properties in D=l+1 for arbitrary SU(N) and present the 
results in D=l+1, SU(3). 



2 Cutoff method 

The cutoff method [13] consists of numerical analysis of the hamiltonian in the occupation 
number representation. First we introduce the bosonic and fermionic creation and annihilation 
operators at % a , a*, f° i.e. 

< = -j=K + ipi), K, a^ b ] = S^6 ab , {/:, /tf } = 8^8 ab \ 
Next we truncate the Hilbert space to the maximal number o quanta 

TlB = a*l a b> n B < nBmax, 

l,b 

compute matrix elements of H and diagonalize the resulting finite matrix. In this way one can 
analyze the spectrum dependence on a cutoff riBmax- There is a dramatic difference between 
the behavior of the continuous and discrete spectrum with cutoff. Namely 

E^ max = E m + 0(e~ nBmax ) - discrete spectrum, 
finBmax _ q( j _ con tinuous spectrum, 

^Bmax 

where m is an index of the energy level m = 1, . . . , nsmax + 1- The limit nsmax — > oo is called 
the continuum limit. In the case of the discrete spectrum case the energy levels converge rapidly 
to the exact eigenvalues of the hamiltonian. This may not be surprising, however it is interesting 

1 There are several choices of fermionic /t" operators. Since we do not make any 
explicit calculations here we refer the reader to [13] for details. 
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to see how fast is the convergence. For details the reader is referred to [14]. In the continuous 
spectrum case things are different. The convergence is very slow and all the eigenvalues vanish 
in the infinite cutoff limit. In the continuum limit the spectrum is continuous and the only way 
to restore it from cut Fock space is to put the following scaling [15] 



m(n Bm ax) = const. y/n Bm ax ^=> E ^f(nB^) ~* E - ( 2 ) 

It was claimed in [15] that this scaling law should work independently of the theory whenever 
one can define scattering states asymptotically. The argument for the above claim is based on the 
following fact. The eigenvalues of the momentum operator in ordinary d=l quantum mechanics 
in cut Fock space are zeros of Hcrmite polynomials H nBmax (x) the asymptotic behavior of which 
is . 1 [14,15]. Therefore, once the momentum operator is defined, its spectrum cutoff 

dependence should be 1 for large riBmax- 

The E^ ma:c values for fixed nsmax give the opportunity to calculate regularized (nsmax 
dependent) Witten index. If the spectrum of the supcrsymmetric hamiltonian H is discrete then 
the index counts the difference between bosonic n° and fermionic n° ground states i.e. 

I w = Tr{-l) F e-P H = YS-V nm) e~ f>Em =<~ n % 

m 

where F is a fermion number. This quantity is (3 independent. The cutoff makes it (3 and n Bmax 
dependent i.e. 

Q 9 (P,n Bmax )= J2 {-l) F ^e-^- Bma \ (3) 

m— 1 

If the spectrum of the hamiltonian H is continuous then the I\y depends on (3 and the 
difference n® — may be obtained by taking the f3 — > oo limit. On the other hand the f3 — > 
limit is easier to compute, therefore one introduces the boundary term 5Iw using the following 
trick [9] 



f°° d 

6I W = I w (oo) - I w (0) = dp—I w {p). 

Jo a P 



2.1 D=l+1,3+1,9+1 SU(2) SYMQM 

In D = 1 + 1 case the hamiltonian H = \p a p a + gx a G a , where G a is the SU(N) generator, is 
free in a gauge invariant sector. There are as many fermion sectors as the grassmann algebra 
allows i.e. 1 boson sector and iV 2 — 1 fermionic sectors. Since the gauge group is SU{2) we will 
denote them as | F = 0), | F = 1), | F = 2), | F = 3). We also have the particle-hole symmetry 
which relates sectors | 0) 3) and | 1) <-»| 2) hance the analysis of the first two sectors is 
sufficient. There is also supersymmetry which relates sectors | 0) <->| 1) and | 2) <->| 3) therefore 
the whole information about the spectrum is in fact in the first sector. Supersymmetry does not 
communicate between sectors | 1) and | 2) which is exceptional for SU(2). Since the particle-hole 
symmetry relates sectors with different fermion number, it is evident that the regularized Witten 
index of this model vanishes. It is however interesting to compute the restricted Witten index 
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which is defined in first two sectors only and the exact answer is \ [16] which was also confirmed 
numerically. 

In D=3+l dimensions the hamiltonian (1) is not free due to the quartic potential term. 
There are 6 fermionic sectors. The particle-hole symmetry relates sectors with the same fermion 
number hance the eigenstates from these sectors do not cancel under the sum (3). The analysis 
of the index [10] shows that in this case 

I w (oo) = I w (0) +61 = -- - = Witten index for D=3+l, SU(2). 

On the contrary the cutoff analysis gives the non zero value [17]. The index converges towards 
i which is exactly the value of the iV(0) not I\y (oo). It seems that the cutoff method somehow 
does not contain the boundary term SI. 

This model is the first non-trivial one where the scaling (2) was confirmed i.e. the spectrum 
of a free particle p 2 /2 can be recovered provided eqn. (2) is applied. Moreover, in fermion 
rich sectors both discrete and continuous spectrum is present which precisely corresponds to 
conclusions of [7]. 

The analysis of the supermultiplets is even more interesting. Each eigenstate is labelled by 
three quantum numbers: energy E, angular momentum 1 and fermion number F. Therefore each 
state can be represented by a dot in R 3 space. It can be proved [17] that supersymmetry links 
these dots in such a way that the emerging geometrical object representing each supermultiplct 
is a diamond. This picture very nicely catalogues all the supermultiplets and it is independent 
of a gauge group. 

In D = 9 + 1 dimensions case we only note the astonishing difficulties that emerge [18]. Since 
we have the SO(9) symmetry the second order Casimir operator is 

J 2 = J ik = x\v k } + \& k ^ a , ^ k = -J[F,r fc ]. 

i<k 

Normally we would have expect the SO (9) singlet to be the Fock vacuum | 0). This is not he 
case here since one can prove that J 2 | 0) = 78 | 0) [18]. The empty state is not invariant under 
rotations! This is a surprising fact and it means that the SO (9) singlet is somewhere else. Where 
is it? The model has 24 fermionic sectors and it was found that the singlet happens to be just 
in the central F=12 sector. 

3 The general properties ot the D=l+1, SU(N) SYMQM 

Since the eigenstates in SYMQM are the gauge singlets therefore it is reasonable to ask about 
the convenient SU(N) invariant basis. It is evident that states belonging to such basis have to 
be of the form 

T hc ... de ..a\a\...flfl... | 0), (4) 

where Tb c ...de... is some SU(N) invariant tensor made out of structure tensors f a bc,d a bc,Sab- We 
now proceed to chose linearly independent states from (4). 
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3.1 Birdtracs 



In order to deal with the variety of all possible tensor contractions we introduce the diagrammatic 
approach (figure 1). 




Figure 1: Diagrammatic notation of invariant tensors. 

Each leg corresponds to one index and summing over any two indices is simply gluing appro- 
priate legs. Structure tensors fij k , dijk are represented by vertices and Sij is a line. Any tensor 
may now be represented by a graph. Such diagrammatic approach has already been introduced 
long time ago by Cvitanovic [19]. In, general one can construct loop tensor which by definition 
is a tensor that diagrammatically looks like a loop however it can be proved [20] that any such 
loop can be expressed in terms of forests i.e. products of tree tensors ( figure 2) 

W 24 4 3 4 3 2 3 3 2 

12 13 12 12 14 14 

Figure 2: An example of loop reduction for a square made out of tensors. 

Therefore we are left with tree tensors only. These however can be easily expressed in terms 
of trace tensors Tr(T a Tb ■ ■ .) where T a are SU(N) generators in fundamental representation. 
With the use of the following matrices = alTb, = /^T(, any gauge invariant state can be 
obtained by acting with an appropriate linear combination of products of trace operators 

Tr{A^F^A^FK..A^F^), 

on Fock vacuum | 0) . Due to the grassmann algebra the number of F matrices under the trace 
cannot be grater then TV 2 — 1 i.e. k < N 2 — 1. Moreover the Cayley-Hamilton theorem for A 
matrices gives ik < N. The remaining set of states is still linearly dependent and the further 
analysis requires separate study of each SU (N) . The basis states in F=0 sector are of the form 

| i 2 ,i 3 *jv> =Tr i2 (A* 2 )Tr i3 (A^)...Tr i »(Ai N ) | 0). 

We see that there are as many states with given number of quanta ns as there are natural 
solutions of the equation 2« 2 + 3h + • • • + Ni^ = Ub- For U(N) this would be exactly p(ub) 
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- the partition number of Ub- For SU(N) this is a little less then p(ns) however it still grows 
exponentially with ns- 

In order to solve the model in bosonic sector one has to compute the following scalar product 

which in principle is a tedious, but not impossible, task . 

Let us discuss, the "bilinear" basis which by definition is the following restricted SU(N) basis 

| 2n) = {A^A^) n | 0), At) = at at (5) 

which was introduced in [2] in D=3+l case. In this basis the non zero hamiltonian matrix 
elements are easy to derive. First we write the commutation relations 

N 2 - 1 

[(AA), At)™] = ^(AtAtf-^AtA) + An{n - 1 + — - — )(A t A t )"- 1 , (6) 

[(AAt), (At At)™] = 2n(AtAt)". (7) 
Using (6) we obtain norms for | 2n) i.e. 



N 2 - 1 

cL : = i 2n I 2n) = 4n(n - H — )c£_ 1; c 2r , 



™ N 2 — 1 
J]; 4fe(fe - 1H ^ — ), c = l. 

\ fe=i 



In the orthonormalized basis | 2n) = ^- | 2n) the non vanishing matrix elements of the hamil- 
tonian 

H = l - PaPa = - J((AtAt) + (AA) - 2 (At A) - (N 2 - 1)), 

are 

N 2 - 1 

(2n \H\ 2n) = n -\ 

(2n + 2 \H\ 2n) = (2n \H\ 2n + 2) = ~ y/ (n + l)(n + N 2 - 1). 

Therefore it is straightforward to proceed with the cutoff analysis (figure 3). We see that thee 
is no quantitative difference between SU(2) and eg. SU(100) case. This is not what we have 
expected and it means that the restricted basis (5) simplifies too much. 



and 
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Figure 3: The cutoff dependence of spectrum for SU{2) and SU (100) in 
"bilinear" basis. 



4 D=l+1, SU(3) SYMQM 

Here we present the calculations of Hamiltonian matrix elements in a complete basis in bosonic 
sector. The basis vectors and the scalar products, we are interested, in are 

| M> = (AUtftAUUt)' | 0), tf i , j jl =(i,j\i',j l ). (8) 
^ At) = d ijk a\a)a\ 

The only non vanishing elements of S 1 ^ are the ones obeying the constraint 2i + 3j = 2i' + 3f. 
Therefore it is convenient to work with the following symbol 

W? j = (i,j\{AAA) 2k {AtAt) 3k \i,j), 

which has the advantage of reproducing all non vanishing TV^-^'s. It is tedious but possible to 
obtain formulas and recurrence equations for W^. We shall omit the lengthy derivation and 
only give the results. 

First we solve the recurrences for IT™ and WK. We have 



< = 96fc(2fc-l)(9fc 2 ~l)(9fc 2 -4)IT fc c T 1 , W* = A(3k + i)(3k + i + 3)Wt 1 , < = 1. (9) 

therefore (9) gives an exact formula for W k . The term is computed from the following 
recurrence 

W Q ) = a jk W^ + & fc W&._ 2 + 7 jk W*+\, 

where 

a jk = 48(2fc + j)(2fc+ j - l)(3fc- l)(3fc- 2)(3fc + 3j + 2)(3fc + 3j + 1), 

jk = 72(2fc + j)(2k + j - l)j(j - l)(9fc 2 + 9kj - 2), 

ljk - 27(2fc + j)(2k + j - l)j(j - l)(j - 2)(j - 3). 

This recurrence stops on given by (9). The general term is now computed from yet 
another recurrence 

W£ - 4(t + 3fe)(t + 3fc + 3j + 3)W?_ij + 3j(j - 1)W££_ 2 , 



which stops on Wqj and W^. The whole norm matrix (8) can now be computed. It should be 
noted that the elements of the N matrix were obtained independently by computing the scalar 
products | i',f) with use of the program written in Mathematica [13]. In this way all the 
recurrences presented here were confirmed up to tib = 12 i.e. for such that 2i + 3j < 12. 
This matrix is in fact the Gram matrix which indicates that we still have to orthogonalize the 
basis. We will not do so however. In order to represent the hamiltonian H in orthogonal basis 
we follow [21]. It is sufficient to calculate the Gram matrix G and proceed with the following 
similarity transformation 



Figure 4: The cutoff dependence of spectrum in D = 1 + 1, 51/(3), F = 0. 

It is clear that the spectrum seems to be far more complicated than in SU (2) case. The lines 
in figure 4 are divided into groups where they converge together. This can be understood in the 
following way. In SU(3) we have two Casimir operators T a T a and d a bcT a TbT c where T a 's are 
SU(3) generators. In cut Fock space the second one does not commute with the hamiltonian 
therefore the cutoff hb breaks the SU(3) symmetry. In — > oo limit the symmetry should be 
restored which corresponds to grouping of the lines in figure 4. 

5 Summary 

SYMQM models reveal variety of application in several areas of physics ( Yang-Mills theories, 
supersymmetry, strings) hance their detailed analysis is of interest. Although they are rich in 
symmetries ( SU(N) , SO(d), supersymmetry ) the exact solutions are missing in the literature 
forcing one to apply numerical methods. The cutoff method presented here is working surpris- 
ingly well, however to get any of results of sections [2,3,4] one had to employ a lot of theoretical 
work which in some cases gave exact results (e.g. the structure of supermultiplets) . The analysis 
of D=l+1 SYMQM for arbitrary SU(N) is very encouraging and gives a hope to proceed with 
the N — > oo limit. 



H ort = G-*HG-t. 

The results of the cutoff analysis are presented in figure 4. 




20 



60 



8 



6 Acknowledgments 



I thank the organizers of the XLV'th Cracow School of Theoretical Physics and the organizers of 
the XXIIIrd International Symposium on Lattice Field Theory for their kind invitation and sup- 
port. I also thank J. Wosiek for discussions. This work was supported by the Polish Committee 
for Scientific Research under grant no. PB 1P03B 02427 (2004-2007 ). 



References 



[i] 

[2] 
[3] 
[4] 

[5] 
[6] 
[7] 
[8] 
[9] 
[10 

[11 
[12 

[is; 

[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 



M. Claudson and M. B. Halpern, Nucl. Phys. 250 (1985) 689. 6189; . 
M. Liischer, Nucl. Phys. B188 (1981) 513. 6189; . 

E. Bergshoef, E. Sezgin, P. K. Townscnd, Ann. of Phys. 185 (1988) 330; 

T. Banks, W. Fischler, S. Shenker and L. Susskind, Phys. Rev. D55, (1997), 6189; 
|hep-th/9610043| 



S. Samuel, Phys. Lett B411, (1997), 268; hep-th/9705167 
B. Simons, Ann. of Phys. 146, (1983), 209 

B. de Wit, M. Liischer and H. Nicolai, Nucl. Phys. B320 (1989) 135. 
M. B. Halpern, C. Schwartz, Int. J. of Modern Physics 13 (1998) 4367; |hep-th/9712133 
P. Yi, Nucl. Phys. B505 (1997), 307; |hep-th/9704098 



S. Sethi, M. Stern, Commun. Math. Phys., 194 (1998) 675, hep-th/9705046 
M. B. Green and M. Gutperle, |hep-th/9711107 

G. Moore, N . Nekrasov, S. Shatashvili, Commun. Math. Phys. 209 (2000) 77, 
hep-th/9803265 

J. Wosiek, Nucl. Phys. B644 (2002) 85-112; |hep"-th/0203116| 

M. Trzetrzelewski, Acta Phys. Polon. B35 (2004) 1615; |hep-th/ 0407059 

M. Trzetrzelewski and J. Wosiek, Acta Phys. Polon. B35 (2004) 1615; |hepTh /0308007. 



M. Campostrini and J. Wosiek, Phys. Lett. B550 (2002) 121-127; hep-th/0209140 
M. Campostrini and J. Wosiek, Nucl. Phys. B703 |hep-th/0407021 



J. Wosiek, Phys. Lett. B619 (2005) 171-176 ; hep-th/0503236 



P. Cvitanovic, Group Theory, http://www.nbi.dk/GroupTheory/ 



M. Trzetrzelewski, J. Math. Phys. 46 , 103512 (2005) ; hep-th/0505084 



P. van Baal, The Witten Index Beyond the Adiabatic Approximation, in: Michael Marinov 
Memorial Volume, Multiple Faces in Quantization and Super symmetry, eds. M. Olshanetsky 
and A. Vainshtein ( World Scientific, Singapore, 2002), hep-th/0 112072 



9 



